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Abstract 

We investigate the exclusive photoproduction of a heavy timehke photon which 
decays into a lepton pair, jp — > p. This can be seen as the analog of deeply 
virtual Compton scattering, and we argue that the two processes are comple- 
mentary for studying generalized parton distributions in the nucleon. In an 
unpolarized experiment the angular distribution of the leptons readily provides 
access to the real part of the Compton amplitude. We estimate the possible size 
of this effect in kinematics where the Compton process should be dominated by 
quark exchange. 
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1 Introduction 



A considerable amount of theoretical and experimental work is currently being devoted to 
the study of generalized parton distributions, whose measurement could make important 
contributions to our understanding of how quarks and gluons assemble themselves to 
hadrons 0, |^. The theoretically simplest and cleanest of the exclusive processes where 
these distributions occur is deeply virtual Compton scattering (DVCS), i.e., 7*p — >• 7p 
in kinematics where the 7* has large spacelike virtuality while the invariant momentum 
transfer t to the proton is small. In the present paper, we investigate the "inverse" process, 
IP ~^ 1*P small t and large timelike virtuality of the final state photon. We shall refer to 
this as timelike Compton scattering (TCS). This reaction shares many features of DVCS, 
although the timelike character of the virtual photon entails some specific differences. 
The combination of data on DVCS and TCS would offer a powerful tool to make sure we 
understand the reaction mechanism, and eventually to obtain stronger constraints on the 
generalized parton distributions than DVCS alone would provide. 




Figure 1: Real photon-proton scattering into a lepton pair and a proton. ^ stands for an 
electron or a muon. 

The physical process where to observe TCS is photoproduction of a heavy lepton pair, 
IP ~^ A* V~ p or ■yp —>■ e"^e~ p, shown in Fig. |I]. Despite the close analogy to real photon 
production ep — ^ e^yp or fip — > fijp, where DVCS can be accessed, the phenomenology of 
these reactions shows important differences. In both Bethe-Heitler (BH) mecha- 

nism contributes at the amplitude level. Contrary to the case of DVCS, this contribution 
always dominates over the one from TCS in the kinematical regime where we want to 
study it. On the other hand, the interference between the TCS and BH processes can 
readily be accessed through the angular distribution of the lepton pair, whereas the cor- 
responding observable for DVCS is the lepton charge asymmetry and requires beams of 
both positive and negative charge. 

This paper is organized as follows. In section |^ we review the kinematics, factorization 
properties, and helicity structure of the Compton amplitude in the general case where the 
two photon virtualities are different, but at least one of them is sufficiently large to provide 
a hard scale. In section ^ we discuss specific features related to the timelike nature of the 
outgoing photon in TCS. We develop the phenomenology of exclusive photoproduction 
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Figure 2: Handbag diagrams for the Compton process (|l]) in the scahng hmit. The 
plus-momentum fractions x, ^, rj refer to the average proton momentum |(p + 



of a lepton pair in section ^, taking into account the Bethe-Heitler and the Compton 
processes and their interference. In section ^ we present estimates of cross sections and of 
asymmetries suitable to extract information on the Compton signal. Section ^ contains 
our conclusions. In an appendix we discuss the relevance of parton densities at very small 
X when modeling generalized parton distributions with a double distribution ansatz. 

2 The Compton amplitude 

Both DVCS and TCS are limiting cases of the general Compton process 

i*{q)+p{p)-^i*{q')+p{p'). (1) 

where the four-momenta q and q' of the photons can have any virtuality. We will also use 
A = p' — p, the invariants 

Q' = -q\ Q'^ = 4\ s = {p + qr, t = A^ (2) 

and write M for the proton mass. In the region where at least one of the virtualities is 
large, the amplitude is given by the convolution of hard scattering coefficients, calculable 
in perturbation theory, and generalized parton distributions, which describe the nonper- 
turbative physics of the process. To leading order in one then has the quark handbag 
diagrams of Fig. 0. The arguments for factorization given in based on the analysis of 
Feynman graphs, hold both for large spacelike and for large timelike virtualities We 
thus define the scaling limit as + ^ oo at fixed t and fixed ratios q'^/s and q''^/s. 
For our subsequent discussion let us recall the expression of the hadronic tensor 

= d'xe-''^-^p{p')\TJ^^{x)jL{0) \p{p)), (3) 

^In contrast, approaches based on the operator production expansion [Q require {q+q'y to be spacehke, 
which is equivalent to + q'^ < when t can be neglected. 
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where eJ^^{x) is the electromagnetic current with e denoting the positron charge. In the 
scahng hmit we have to leading order in 



9t 



u{p). 



(4) 

This expression holds in reference frames where both proton momenta p and p' have 
small transverse components of order and are moving fast to the right, i.e., have 
large plus-components. Light-cone coordinates are defined as = (f° ± v^) /V2 for 
any four- vector v. The transverse tensors Qt and have as only nonzero components 
—g^ = —g^ = = —e^ = 1. Following the notation of we have introduced the 
convolutions 
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^ — X — ie ^ + X — ie' 

t, — X — te t, + X — le'^ 

of the generalized quark distributions defined in summed over quarks of flavor q and 
electric charge ecq. The scaling variables ^ and t] are given by 



2{p + p') ■ {q + q') 2s + g2_g/2' 



(g - q') ■ (g + q') ^ Q' + Q" . . 

{p + p') ■ {q + q') ^ 2s + Q2_Q'2' I ) 

where the approximations hold in the kinematical limit we are working in. x, ^, and r] 
represent plus-momentum fractions 

{k + kT {q + q'Y (p-pT 

{p + pT' {p + pT' '^{p + p'r' ^' 

The expressions (|D and (|^) reveal that the two-photon amplitude is independent of the 
photon virtualities at fixed ^, t] and t. In the case of spacelike q = q' this is just Bjorken 
scaling. To be precise, the independence on and g'^ only holds up to logarithmic 
corrections: the photon virtualities provide the hard scale of the process and thus enter 
through the factorization scale dependence of the parton distributions, which we have not 
displayed above. The corresponding evolution equations are well-known |^, |, 0, and 
as usual we will refer to — 1 < x < —t] and rj < x < 1 as the DGLAP regions, and to 
—rj < X < 1] as the ERBL region of the parton distributions. 
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Let us now recall the helicity structure of the two-photon process in the scaling limit. 
Contracting the hadronic tensor with polarization vectors e of the incoming and e' of the 
outgoing photon, one obtains the helicity amplitudes of (|^) as 



where A (A') denotes the helicity of the incoming (outgoing) proton and (/x') the 
helicity of the incoming (outgoing) photon. Parity invariance provides the relations 
J^-\'-^'-x-^, ^ (•_i)A'-m'-a+m j^//AV,AM_ From (g) one easily finds that the quark hand- 
bag diagrams only generate helicity conserving transitions between transverse photons, 
jy-A'+,A+ M^'-'^-. At order as one further has amplitudes M^'°'''^°, provided of course 
that both photons are off-shell §]. Double helicity flip amplitudes and 
are generated at order by gluon transversity distributions Finally, transitions 

involving one transverse and one longitudinal photon are suppressed by one power of 
the large scale Q or Q'. These twist-three contributions^ have been studied in [jll 



and 



These studies were 



twist-four contributions to the double helicity flip amplitudes in |]T2 
performed for large spacelike virtualities; whether they can be extended to the timelike 
case is a question beyond the scope of this paper. 

The DVCS and TCS processes are limiting cases of (|I]) where one of the photons is on 
shell. From (|^) we readily see that to leading-twist accuracy one has C, = rj in DVCS and 
^ = —r] in TCS. The convolutions (^ obey 



ni{-r],r],t) = [ni{r],T],t) 
£i{-r],r],t) = \^£i{ri,r],t) 

which leads to the simple relations 



ni{-r],T],t) = - [ni{7],r],t) 



A' + ,A+ 
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TCS 







DVCS' 
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TCS 



A' + ,A+ 



DVCS 



(9) 



(10) 



between the helicity amplitudes for TCS and DVCS at equal values of r] and t. These 
relations should be evaluated at corresponding values of Q'^ and since the photon 
virtualities play analogous roles in providing the hard scale of the respective processes 
and thus enter in the scale dependence of the parton distributions. The relations (p!OD 
tell us that at Born level and to leading twist one obtains the amplitudes for TCS from 
those of DVCS by changing the sign of the imaginary part and reversing the photon 
polarizations. To this accuracy, the two processes thus carry exactly the same information 
on the generalized quark distributions. 

We remark that the relations (§) and hence (|10D no longer hold at 0{as), neither for 
the one-loop corrections to the quark handbag diagrams in Fig. ^ nor for the diagrams 
involving gluon distributions. On general grounds, the phase structure of the two processes 
is in fact different. Whereas the only discontinuity of the two-photon amplitude in DVCS 



^We use here the dynamical definition of twist, where twist n contributions to the Compton amphtude 
are suppressed by n — 2 inverse powers of the large scale. 
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(a) (b) 



Figure 3: The loop momentum configurations x = rj where the Born level amplitude 
receives its imaginary part in (a) DVCS and (b) TCS. Short vertical lines indicate on- 
shell quark lines in the hard scattering, plus-momentum fractions ±2?; and refer to the 
average proton momentum + The corresponding configurations for x = —1] are 
obtained by reversing the charge flow of the quark line. 



kinematics is in the s-channel, the TCS amplitude has discontinuities in both s and Q'^, 
with one-loop hard scattering diagrams contributing to both cuts. In situations where 
0{as) contributions are important, the DVCS and TCS processes will have a different 
dependence on the generalized parton distributions. TCS and DVCS together can then 
constrain them more effectively than either process alone. The detailed study of TCS at 
one-loop level is beyond the scope of this work, and we will base our numerical studies on 
the Born level expression (H). 

It is worthwhile to compare the momentum configurations in DVCS and TCS from 
which the Born level convolutions receive their imaginary parts. From Fig. ^ we see 
that in both cases there is a quark line with zero plus-momentum coming from the proton, 
and that in both cases it is attached to the real photon, i.e., to the final state in DVCS 
and to the initial state in TCS. 

We conclude this section by defining the variable 

Q" Q" .... 

^ 2p-q s-M^ ^ ' 

for the TCS process as the analog of the Bjorken variable xb = Q"^/ (2p • q) in DVCS. The 
similar roles played by these quantities reveals itself in their relations with t], which to 
leading-twist accuracy reads r/ = r/(2 — r) for TCS and r] = xb/{'2 — xb) for DVCS. 

3 The timelike photon 

Processes involving timelike photons can have markedly different features than processes 
controlled by large spacelike virtualities. These features usually do not arise to leading 
order in perturbation theory, which is the approximation we will work in here. A closer 
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(a) (b) 

Figure 4: Bom level diagrams for the (a) Drell-Yan process pp X and (b) e~^e 

annihilation into hadrons, e~^e~ X. 



look at the Born level diagrams reveals nevertheless important similarities and differences 
between timelike processes, which we now briefly discuss. 

The reaction which at first sight is most similar to TCS is Drell-Yan pair production 
in hadron-hadron collisions. In that case, the 0(0;^) corrections to the Born graph of 
Fig. ^ have considerable size and make up for most of the much discussed /^-factor of 
this process. A way to understand them is the occurrence of large contributions enhanced 
by TT^, which can be traced back to the correction of the quark-photon vertex for spacelike 
7* and on-shell quarks [13|. Notice that in the TCS Born graphs of Fig. ^ only one of 



the two quark lines attached to the 7* is on-shell, whereas the other one is off-shell by 
order Q'^. One might argue that the second line does become on-shell in the imaginary 
part of the amplitude, as indicated in Fig. ^b, but there is an important difference: the 
quark lines in the Drell-Yan diagram and one of the lines in TCS physically correspond to 
small virtualities as they are directly attached to parton distributions, i.e., to quantities 
describing long-distance physics. This is not the case for the vertical quark line in the 
TCS diagrams. Technically, the singularity of its propagator can be avoided by analytical 
continuation of the loop momenta, whereas the singularities associated with the lines 
attached to a parton distribution are pinched 0, The analogy between the two 

processes must hence be used with care, and in particular one cannot easily infer on the 
size of the 0{as) corrections from the experience with the Drell-Yan process. 

A second issue in processes with timelike photons is the importance of resonance 
effects, which are beyond the realm of perturbation theory. At invariant photon masses 
above 4 or 5 GeV, excluding of course the region of the T resonances, the comparison of 
leading-twist perturbative calculations and data works rather satisfactorily for the Drell- 
Yan process, cf. the data compilation in |T^. The situation for masses below the J/"^ is 
difficult to assess, mainly due to background lepton pairs from the weak decays of b and c 
quarks |T^. This type of background does of course not affect TCS, where we are dealing 
with exclusive lepton pair production. As for inclusive e~^e~ annihilation into hadrons, the 
recent BES data in the mass region from 2 to 3 GeV is remarkably fiat and close to 



the leading-twist result. The same holds for the data above 5 GeV, cf. e.g. [|18[, excluding 



again the T region. Between 3 and 5 GeV on the other hand, resonance structures are 
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clearly visible [ff7 



Again one should keep in mind that the importance of resonance effects may be differ- 
ent in all these processes. In line with our above analysis, we remark that in the tree level 
diagram for inclusive e^e~ annihilation, Fig. ^o, both quark and antiquark correspond 
to large virtualities. Technically, the cross section is calculated as the imaginary part 
of the photon vacuum polarization, where the quarks appear in a loop and are indeed 
far off-shell. We notice that in both Drell-Yan production and e^e~ annihilation one 
has quark-antiquark configurations with comparable virtualities. In contrast, we have 
asymmetric configurations in TCS, with one quark line soft and the other far off-shell. 
Furthermore, the space-time structure of TCS is such that the 7* is formed from a gg-pair 
only in the ERBL region of the parton distributions, while in the DGLAP region the 
parton- level process is photon radiation off a quark or antiquark, q —>■ '~f*q ot q '~f*q. 

To conclude, we estimate based on e'^e" —>■ X and the Drell-Yan data that ranges of 
Q' where the leading-twist description of TCS may work should be between about 1.5 to 
2 GeV and the J/\E' mass, and above the charmonium resonances. We stress however that 
the reaction mechanism in the TCS process displays important differences, and that one 
will have to see in the data how parton-hadron duality manifests itself here. 



4 Observing TCS in lepton pair production 

4.1 Some kinematics 

Let us now specify the variables we use to describe the lepton pair production process 
depicted in Fig. |I], in addition to those already introduced at the beginning of Sect. ^. 
A useful quantity is the transverse component of the momentum transfer A with 
respect to p and q in the •yp cm. It is related to the scattering angle 0cm in that frame 

by 

sme,^ = —^, (12) 



where A^ = |At| and r = ^ {s - Q''^ - M'^f - AQ''^M'^. In the limit of large Q'^, large 
s, and small — t, we then have 

-t.^--^ (13) 
up to relative corrections of order M^/Q'^. For the lepton pair, we use the lepton velocity 



/? = - 4m2/g'2 (14) 

in the cm., where mi denotes the lepton mass. In the same frame we introduce 

the polar and azimuthal angles 9 and 99 of fc, with reference to a coordinate system with 
3-axis along —p' and 1- and 2-axes such that p lies in the 1-3 plane and has a positive 
l-component.[] This is shown in Fig. |^. In terms of Lorentz invariants, our angles are 

^They correspond to the decay angles 9 and for vector meson photoproduction introduced by 
Schilling, Seyboth and Wolf with their vector 7? along k. 
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Figure 5: Sketch of the kinematical variables and coordinate axes in the 7p and cm. 
frames. Notice that the coordinate systems differ from the one we used in the Compton 
amphtude (p, where p and p' have positive 3-components. 



given by 



2{k-k')-p 
2{k-k') ■ [p-p') 



Prcos9, 
a(3 



\ 



2{s - M2) Q'At 



cos 6 



— p sm U cos (p, 



4e^''^" pJukX = (^{s- M^) Q'At sin 6 sin if, 



(15) 



where our convention for the completely antisymmetric tensor is 69123 = 1, and the sign 
factor 0" = ±1 is determined by 



\ 



(g'2 - ty - 



2{s - M^)Q'At 



I2\ 



(16) 



The form of the second equation in ([T5|) is useful in our kinematics, where At is small 
and cr = 1. 

As polarization vectors e(A) for the incoming photon we take e(±) = {Te^^'^ -ie^'^^) / \/2, 
where e*-^-* and e^^-* are unit vectors along the 1- and 2-directions in the 7p cm. as shown 



in Fig. H. Our polarizations e'(A') of the outgoing photon are e'(±) = (=Fe 



'(i)_,e'(2) 



and e'(0) = e''-'^^ with unit vectors along the coordinate axes in the cm. described 
above. 



4.2 The Bethe-Heitler contribution 

The Bethe-Heitler amplitude is readily calculated from the two Feynman diagrams in 
Fig. ^ We parameterize the photon-proton vertex in terms of the usual Dirac and Pauli 



9 



Figure 6: The Feynman diagrams for the Bethe-Heitler amphtude. 



form factors Fi{t) and F2{t), normahzing F2{0) to be the anomalous magnetic moment of 
the target. We find for the BH contribution to the unpolarized jp cross section 



da 



BH 



/5 



(^1 



t 



4M2' 



A 



B 



(s - M^)^A^ -ta{a + h)- - t {AM' - t)Q 
+ ^ 
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dQ''^dtd{cose)dip 47r(s - M2)2 -tL 
where we have used the abbreviations 

{ -t){a + b)-is-M^)bY + t (4M2 - t) - t) 



(17) 



B 



iQ" + ty + b' + 8miQ 



The cross section depends on the angles 9 and through the scalar products 
a = 2{k-k')-p', b = 2{k-k') -{p-p') 



(19) 



given in Eq. (15) above, and through the product of the lepton propagators in the two 
BH diagrams, 



{q - kf 



{q-k 



l\2 



iQ'^ - tf - b^ 
4 



(20) 



These expressions are rather lengthy, but simplify considerably in kinematics where t, 
and m| can be neglected compared to terms going with s or Q'^. We then have 



and 



da 



BH 



at 



1 l + cos^^ 



sin" 9 



dQ'^ dt d{cos 9) dip l-ns^ -t sin^ 9 



(^1 



4M2' 



(21) 



• (22) 



We see that the product L of lepton propagators goes to zero at sin 6' = in this approx- 
imation. Closer inspection reveals that when sin 6' becomes of order At/Q' or m^/Q' the 
approximations (pTf) and (E^) break down and one must use the full expressions. 
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Let us see how small the product L can become. At fixed s, Q'^, t, we find with 
Eqs. (IT5|) and (EDI) that L assumes a minimum value 



,2 



for 



L^r. ^ Q"m] + Q'^Al (23) 

(1 -r)^ 

2 At- cos 09 , 

tan^^i,^ -— f (24) 
Q 1 — r 

up to corrections of order t/Q'"^, M^/Q'^, m\IQ''^ . For ~ the leptons t~ and are 
nearly coUinear with the initial photon in the 7p cm. They have transverse momenta of 
order Ay with respect to p and q and share their total longitudinal momentum in a highly 
asymmetric way. In our numerical studies we will impose a cut on Q which ensures that 
L remains of order Q'^, thus staying away from the region where the BH cross section 
becomes extremely large. 

We finally remark that as long as L is of order Q"^ the terms going with 1/L in ([TsD 
are suppressed at least like rn^Q'^/L compared with the leading behavior of A and B. 
For a large range in 9 the BH cross section (|l^ will thus approximately behave like 1/L 
instead of 



4.3 The Compton scattering contribution 

We now investigate the TCS contribution to lepton pair production. In order to under- 
stand the basics of its interplay with the BH process it is sufficient to consider the leading 
behavior of the Compton amplitude in 1/Q' and in a^, which we discussed in Sect. H. We 
will thus in particular discard 7p amplitudes that change the photon helicity. In line with 
neglecting power suppressed effects in the Compton subprocess we will also drop mass 
corrections of order M^/Q'^ and m|/Q'^ in kinematics and phase space. To this accuracy, 
the contribution of TCS to the unpolarized cross section of 7p t^l' p reads 

dares _ 1 1 + cos^ ^ y_ . 3 



dQ'^dtd{cose)dip 87rs2 g'2 4 



^|^A'-,A-|2_ ^25) 



A,A' 



We note that the (p independence here is a consequence of having neglected photon helicity 
changing transitions. From (^ we obtain 



2a,a 



where TCi, 7ii, Si, £1 are to be evaluated at = rj. Together with Eq. (l2^ ) we see that 
compared with the TCS cross section, the BH contribution is parametrically enhanced by 
a factor Q'^/(— t) and has an extra factor of 1/ sin^6' in the angular dependence. 
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Let us compare the TCS result p5| ) with the contribution of DVCS to the electropro- 
duction process 

m + Pip) ^ ^ik') + 7(g') + Pip'), (27) 

where we have indicated four-momenta in parentheses. Retaining only the leading part 
in 1/Q and of the Compton amplitude, and dropping again mass corrections of order 
M'^/Q'^ and mj/Q"^, we have for the unpolarized cross section 

dapvcs _ «L 1 l + (l-y)' ^,,.y+.A+|2 .281 

Here Sep = (p + kY is the total cm. energy of the ep colhsion, y = iq ■ p)/ik ■ p) the usual 
inelasticity parameter, and if the azimuthal angle between lepton and hadron planes as 
defined in 0. With the relation (p!0D we readily see that to leading twist and leading order 
in as the sums over squared helicity amplitudes in (^) and (^) give identical results for 
corresponding values of r/ = r/(2 — r) and Q'^ in TCS, and r] = 0:^/(2 — xb) and 
in DVCS. To this accuracy, the Compton scattering contributions to the respective cross 
sections only differ by the global kinematic factors given in ( ^51) and (|28|) . 

Comparison of these factors reveals the correspondence between the variables 6 in TCS 
and y in DVCS, which by expressing them in terms of scalar products is found to be 

l + cos6' k ■ p' k-p 1 

(29) 



{k + k') -p' {k- k') ■ p 



where in the first relation we have again neglected mass corrections. At this point we find 
a crucial difference in the phenomenology of the two processes. As is well known p, |20[ the 
relative weight of DVCS and BH crucially depends on y, given that at amplitude level the 
DVCS contribution comes with a factor 1/y relative to the BH contribution. In the region 
of and t defining the DVCS regime, BH dominates for moderate values of y, whereas 
for sufficiently small values of y the Compton contribution wins. Since the quantity 
corresponding to l/y in Eq. (^) is always between —1 and 1, no such enhancement takes 
place for TCS, and we will indeed find numerically that here the BH contribution to the 
cross section is always dominant. The strategy is then the same as in DVCS at moderate 
values of namely to gain information on the Compton process through its interference 
with BH, which can be extracted using symmetry properties of the process. 

Another noteworthy difference concerns the variables r and xb-, which determine the 
values T] where the generalized parton distributions are probed in the two processes. In 
DVCS at fixed collision energy ^Js^p the variables xb and y are not independent since 

= yxB (sep — mf — M^). If at given one wants to vary xb and thus choose a 
kinematical point where to probe the Compton subprocess, one must vary y. In TCS 
on the other hand, one has the relation Q'^ = t {s — M^), independent of the value 
of ^. In order to vary r at given Q'^, one here needs to change the 7p collision energy 
^/s. A continuous spectrum in s is of course automatically obtained if the initial photon 
originates from bremsstrahlung off a lepton beam. 
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4.4 The interference term 



Let us now explore how information on the Compton process can be obtained from the 
interference between the TCS and BH amphtudes. The general strategy is the same as 
described in for the case of DVCS, but we will again encounter important differences 
in the phenomenology of these reactions. 

A key point is that the amplitudes for the Compton and Bethe-Heitler processes 
transform with opposite signs under reversal of the lepton charge. As a consequence the 
interference term between TCS and BH is odd under exchange of the and ^~ momenta, 
whereas the individual contributions of the two processes are even. Any observable that 
changes sign under k ^ k' will hence project out the interference term, eliminating in par- 
ticular the large BH contribution. Clean information on the interference term is therefore 
contained in the angular distribution of the lepton pair. The corresponding observable 
in the electroproduction process ip — ^ ijp is the lepton beam charge asymmetry, whose 
measurement presents important experimental challenges. 

Let us take a closer look at the interference part of the cross section for 7p p 
with unpolarized protons and photons. It is given by 



da 



INT 



1 M 



1 



dtd{cose)dip 



Airs^ -t Q' rVT 



T L 



1 + cos^ e ^ ~ 

COS if — Ke M 

sin 9 



- cos 2^ V2 cos 9 Re M°- + cos 3^ sin 9 Re M+" + O (-^) 



with L and Lq from Eqs. (^0|) and (|2T|) . Here 



~ , At 
ikf"^ = — 
M 



T 



M- 



+ 



At 



F1 + -F2 
. 2 



+ 



r\F, + F2) + 



A. 



T 



2M2" 



A2 

2M2 ^ 



(30) 



(31) 



is the same combination of Compton helicity amplitudes as defined in 0.0 The close 
analogy between TCS and DVCS is manifest, and we see that a 7* with negative helicity 
in TCS corresponds to a 7* with positive helicity in DVCS as we already found in the 
relations (|10D . 

The terms indicated by 0(1/(5') in Eq. (|30| ) have kinematical coefficients suppressed 
by at least one power of 1/Q' relative to the other terms in brackets. Notice that we 
have not approximated the product L of lepton propagators from the BH process. In 
the limit of large Q'^ the factor Lq/ L tends to 1, but we have seen in Sect. ^]2| that this 



approximation becomes increasingly bad as 9 approaches or vr, so that it is useful to 
keep Lq/L in an analysis. The same is true for the lepton propagators in the interference 
of DVCS and BH, as has been emphasized in 0. 



^In contrast to P] our notation here is to list the hehcities of outgoing particles first. With our phase 



convention the transverse polarization vectors of the two photons coincide for = 0, cf. Sect. 4.1. In [p| 



we made a different choice, and the Compton helicity amplitudes here and there differ by an overall sign. 
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We see that without polarization one probes the real parts of the Compton helicity 
amplitudes. Access to the imaginary parts can be obtained with polarized photon beams. 
If the photons have a circular polarization i/, as is the case for a bremsstrahlung beam 
emitted from longitudinally polarized leptons, one has 

da INT da I NT 



dQ'^dtd{cose)d^ dQ'^ dt d{cos 6) dip 



— V 



eq. (30) 

em ^'^ -^0 



aL 1 M 1 L 



47rs2 -t Q' tVT^ L 



1 + cos^ e ^ ~ 

sm (p — Im M 

sin 6 



- sin 2lp V2 cos 6 Im M°" + sin 3^ sin 6 Im + (t^) 



(32) 



The photon polarization dependent and independent terms are simply related by exchang- 
ing sin ^ cos and Im Re . This is not quite the same as for lepton beam polarization 
in the interference between DVCS and BH, where different kinematical factors occur in 
the polarization dependent and independent parts, and where notably the term with 
sin 3^5 Im M"*"" is absent. 

The various terms in the (p dependence of the interference term can for instance be 
projected out by weighting the differential cross section with appropriate functions. The 
weights (L/Lq) cos{n(p) and (L/Lq) sm{n(p) for instance give the terms with cos(n<y9) and 
sm{rnp) in Eq. ( ^O]) and (|32|) , respectively. Notice that these weights are odd under the 
exchange of k and k' and hence do not pick up the BH and TCS contributions to the cross 
section, as discussed above. 

In this way we can project out the various helicity combinations ^ of Compton 
amplitudes, up to relative corrections in IjQ'. Along the lines of this can be used to 
test whether the power behavior in Q' at fixed r and t follows the predictions discussed 
in Sect. |^, i.e., whether arguments based on the large Q'^ limit apply at the finite Q'^ of 
a measurement. If one is in the scaling regime, one can then analyze the photon helicity 
conserving amplitudes in terms of generalized parton distributions. The quark handbag 
diagrams of Fig. |^ give 



(33) 



M l + T] 

where —to = 477^M^/(1 — 1]"^) is the minimal value of —t at given t], up to corrections 
in l/Q'2. 

The above extraction of the the Compton amplitudes requires measurement of the an- 
gle ip. If one integrates the interference term over ip, the photon polarization dependent 
part in (^2]) vanishes because of parity invariance. The integral of the unpolarized contri- 
bution (^) is nonzero, due to the ip dependence of Lq/L and to the (p independent part 
of the terms denoted by 0{1/Q'). This integral can in principle be projected out from the 
cross section because it is odd under 9 ^ n — 9, whereas the BH and TCS contributions 
are even when integrated over cp. The interference signal so obtained is however an order 
1/Q' smaller than what can be seen in the ip dependence of the cross section, and will 
thus be harder to extract. 
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5 Numerical estimates 



In this section we model the generahzed parton distributions (GPDs) and give estimates 
for various observables. We restrict ourselves to moderate values of r and use the leading- 
order handbag approximation (|), (j^) of the Compton amplitude. We omit all terms 
proportional to E"^ and E'^. In the region 0.1 < r < 0.36 and \t\ < 0.4 GeV^ we will 
consider in our estimates, Si is multiplied by kinematical coefficients at most 0.15 times 
those of Hi in (OO) and (|5BD and thus would not significantly change our results. In any 



case, it is at present fairly unclear how to model the distributions E"^, so that taking them 
into account would not improve the reliability of the estimates. As for Si, it is multiplied 
by a tiny coefficient in (p6D and absent in the interference term (SSl). 



5.1 Modeling the parton distributions 



Now we define the model we use for H'^ and if. Following |^ we take a factorizing 
ansatz for the t dependence, 

H^ni^,r^,t) = /i"(x,r/)|Fi"(t), 

KDi^,V,t) = h''{x,v)F''{t), (34) 

with 

Fr(t) = 2Ff(t)+Ff(t), 
Ftit) = Ff(t) + 2Ff(t), 
= F'{t) = gA{t)/gA{^). (35) 

Ff and are the electromagnetic Dirac form factors of the proton and neutron, for which 
we take the usual dipole parameterization p2| . For the axial form factor of the proton we 
take QAit) = QAiO) (1 - t/Ml)-^ with c/^(0) = 1.26 and Ma = 1.06 GeV from [0. For 
strange quarks we make the ansatz 

H'nDix,r],t) = h\x,v)FD{t), (36) 

where the dipole form factor Fz)(t) = {l-t/M^)'^ with My = 0.84 GeV is the same that 
enters in the parameterization of Ff and Fi. Note that via the sum rule for / dx H^{x, 1], t) 
a factorizing ansatz like (|36| ) corresponds to setting the strange quark contribution Fi{t) 
to the Dirac form factor to zero. We remark that several studies [^, |2^ indicate that 
GPDs do not factorize in the simple manner of (0) and (|5B|). The ansatz has however 
the virtue of simplicity and should be good enough for our estimates, as long as we do 
not study the interplay of the r] and t dependence of the cross section. For h'^ and h"^ we 
make an ansatz based on double distributions p6[] . 
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1 rl-x' 

dx' / dy' 

J-l+x' 



6{x — x' — rjy') q{x') 



6{x + x' — rjy') q{x') 



7i{x',y'), 



h'^{x, Tj) 



7i{x',y') 



1 rl-x' 

dx' / dy' 6{x — x' — rjy') Aqv{x') 7r(x', y'), 

J~l+x' 

'\2 „,/2 



3 (1- x'Y - y '' 
4 (1 



(37) 
(38) 
(39) 



We evaluate (|37D with the LO GRV 94 parameterization [^] of the unpolarized distri- 
butions q{x) and q{x), and (|38D with set A of the LO polarized valence distributions 
Aqy{x) by Gehrmann and Stirling ||28|. In both cases we take the factorization scale as 
jj? = b GeV^. We neglect the polarized quark sea, which presently is not well constrained 
by data, and thereby also drop H'^ . In the appendix we shall give a detailed discussion 
of the role played by very small values of x' in the integrals of (^^ and (|38|) , and thus of 
the uncertainties in evaluating them with parton densities only known above some finite 
value of x' . 

Let us stress that the available models of GPDs are fraught with uncertainties, in 
particular in the ERBL region. There, GPDs describe the emission of a qq pair from the 
target, and an ansatz only using the information from usual parton densities should be 
used with care. Dynamical calculations [^, ^ lead in fact to much richer structure in 
the ERBL region than is generated from (^) to (|39|) . Notice also that, while for x > rj 
GPDs are bounded from above [^], no analogous constraints are known in the ERBL 
region. 

A particular type of contribution in the ERBL region is the Polyakov- Weiss D-term 
31 1, which following |32] we take as a flavor SU{3) singlet 



H^x, T], t) = HUx, V, t) = HUx, V, t) = 0(^' 



X 



-D 



X 



(40) 



where G denotes the step function. We make again a factorizing ansatz for the t de- 
pendence, taking the same dipole form factor as in (EBf). For the function D we use the 



parameterization given in equations (23) and (24) of [^], which was obtained by a fit 
to the result obtained in the chiral soliton model p9|. That parameterization is given 
for a factorization scale /i = 0.6 GeV, and we use the leading-order evolution equations 
to evolve it up. Because of mixing we then need the D-term in the gluon GPD of the 
proton, which we take as zero at = 0.6 GeV. Following ^ we take A^^) = 232 MeV 
and A^^) = 200 MeV for the scale parameter in as, switching from 3 to 4 flavors at 
/X = 1.5 GeV. For /i^ = 5 GeV^ we then get 



Diz) 



2.9 Cl^\z) + 0.6 G 



3/2/ 



+ 0.2C5/^( 



(41) 
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Figure 7: The contributions from u, d, and s quarks to ReTii (left) and IniTii (right). 
They are calculated with H'^ = Hl)^ and respectively divided by ^F^(t), Ff(t), and Fu^t). 
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Figure 8: The contributions from u and d quarks to ReTii (left) and ImTii (right). They 
are calculated with our model for H'^ and divided by F'^{t). 



with Gegenbauer polynomials C'^^'^{z). Below, we will give estimates with and without the 
D-term contribution according to ( ^0|) and (^Tj) in order to explore the model dependence 



of our results. 

In Fig. 1^ we show the real and imaginary parts of the convolution integral Hi (—77, r], t), 
calculated from Hod- Decomposing Hi = Hi + Hf + HI we plot the contributions 
from u, d, and s quarks separately. We further divide by appropriate factors |F"(t), 
Ff(t), and Fj^it), so that with the factorizing ansatz (pi]), ( |5^ ) the resulting curves are 
independent of Analogous plots for Hi = Hi + Hf are given in Fig. |[ We observe 
that for r ~ 0.1 the s quark contribution to ReT^i is by no means small compared with 
u and d quarks, although it is tiny in ImHi. This illustrates that at least the real part of 

^Formally, these curves correspond to TY"^'''" at t = 0, which for r 7^ is however outside the physical 
region according to (O). 
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X 



Figure 9: The combination H]^^(x,ri, t) - H^d{-x, r], t) divided by \F^{t), and the com- 
bination H^{x,ri,t) — H^{—x,ri,t) divided by -Fz)(t). Both curves are for r] = 0.11, 
corresponding to r = 0.2. 



the Compton amphtude is not related in a straightforward manner with the usual parton 
densities at x ~ r, given that s{x) only becomes comparable to u{x) and d{x) for x 
considerably below 0.1. 

We do not show in Fig. |^ the contributions from Ho- They are only nonzero in 
the real part, summed over all flavors they amount to a r independent contribution of 
—3.3 Foit) in Tii. The remarkable fact that the D-term contribution to the TCS ampli- 
tude is independent of 1] at fixed t remains true to all orders in perturbation theory. This 
is because due to general scaling properties the hard scattering kernel can be written as 
the leading-order one in (^) times a function of x/t]. Comparing with Fig. ^ we see that 
the D-teim has an appreciable impact on the value of Re Tii in our model. This is sur- 
prising if one compares the functions and if^D themselves. We show this for u quarks 
in Fig. ^, plotting only the charge conjugation even combination H{x,ri,t) — H{~x,ri,t) 
that enters in Compton scattering. One can understand the strong amplification of a 
moderate change in the ERBL region of a GPD by observing that the real part in the 
convolutions (^) is a principal value integral, which involves large cancellations between 
the contributions from |a;| < rj and |x| > 7]. Here is one of the reasons why measuring the 
real part of the Compton amplitude, in DVCS or in TCS, can provide unique information 
on generalized parton distributions. 



5.2 Cross section and angular distributions 



To calculate the TCS amplitude we start with the hadronic tensor (^), evaluated in the 
7P cm. with the 3-axis in the direction of p. In order to preserve gauge invariance beyond 



the leading- twist approximation, we use the prescription of |21] and take 



T 



eq. (4) 



9l 



(42) 
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Figure 10: The BH (solid line) and TCS (dotted line) cross sections for ^/s = 5 GeV and 
Q'2 ^ 5GeV^ integrated over (/? e [0, 27r] and 9 e [7r/4, 37r/4]. The TCS contribution is 
calculated using and H1)d from 



where the index (3 refers to the virtual photon. The subtraction term is formally sup- 
pressed by 1/Q' and has effects of a few percent on the results we will present. The 7p 
cross section is then calculated from (^) and the exact expression of the BH amplitude. 
We have compared the interference term thus obtained with the approximate expressions 
in Sect. 0. For Q'^ = 5 GeV^, \t\ = 0.2 GeV^, and = 5 GeV we find that the approx- 
imation (^) with (P^l ) deviates by at most 10% from what we obtain with (J^). As one 
expects, the situation gets worse for larger values of |t| but improves quickly for larger 
values of Q'^. 

We are now ready to estimate the different contributions to the cross section. In Fig. |1^ 
we show the result for the TCS and the BH contributions to the yj-integrated cross section 
at ^/s = 5 GeV and Q'^ = 5 GeV^. Here and in the following we integrate over 9 from tt/A 
to Sir /A, avoiding the region where the BH contribution becomes hopelessly large. As we 
anticipated in Sect. [4.3| the BH process nevertheless dominates the cross section, with TCS 
contributing less than 5% in this kinematics. In Fig. |10] there is no contribution from the 
interference between TCS and BH since the angular integration selects charge conjugation 
even quantities. Let us therefore investigate the manifestation of the interference term in 
the angular distribution. We restrict ourselves to unpolarized photons here. 

In Fig. 0a we show the dependence of the cross section integrated over 9 in the 
range [7r/4, 37r/4]. With the integration limits symmetric about 9 = 7r/2 the interference 
term is odd under (p n + ip due to charge conjugation, whereas the TCS and BH cross 
sections are even. We separately show the contribution from BH and the sum of BH 
and the interference term. The TCS cross section is fiat in ip to leading-twist accuracy, 
cf. (pSD, and only tiny oscillations are induced by the prescription (|4^) . In the kinematics 
of the figure we get dares / (dQ'"^ dt dip) ~ 0.2 pbGeV"^ when applying the same cut in 9 
and taking the double distribution ansatz (|3^) for the GPDs. 
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Figure 11: (a) The cross section integrated over 6 G [7r/4,37r/4] as a function of (f for 
^/s = 5 GeV, Q'^ = 5 GeV^, \t\ = 0.2 GeV^. The curves represent the BH contribution 
(sohd hne) and the sum of BH and the interference term (dash-dotted hne), calculated 
using Hfj^ and HIjd- (b) The same as in (a) but with the cross section weighted by L/Lq 
before integrating over 6. 



Fig. [TT| b shows the corresponding contributions to the weighted cross section 



dS 



dQ'"^ dt dip 



L(g^y) da 

n/4 Lo{e) dQ'^dtddd^' 



(43) 



We see that the signal is more easily visible after this weighting. The interference term 
behaves now like cosy? up to 1/Q' suppressed terms that are numerically small. The 
weighted BH cross section is almost flat with our cut on 6', in line with our discussion at 



the end of Sect. The TCS contribution is again small here and will not much change 
the picture. As discussed in Sect. ^.4|, information on the interference can in principle also 



be obtained from the yj-integrated cross section. With the same kinematics as in Fig. ^ 
we find that the interference generates an asymmetry in da / {dQ'"^ dt dO) about 9 = 71/2 
which is barely at the 1% level. 

To extract information on the Compton amplitude in a compact way we introduce 



R 



I-2TT 

2 d(p COS(y9 

Jo 



dS 



dQ'"^ dt dip 



2- dS 
dQ''^ dt dip 



(44) 



which projects out the ratio ai/ao of Fourier coefficients in the weighted cross section 
dS / {dQ'"^ dt dip) = J2'^=o CLnCos^rup). Up to 1/Q' suppressed contaminations the numera- 
tor in R is proportional to the combination M of Compton amplitudes , whereas the 
denominator is in our kinematics dominated by the BH part of the cross section. To 
explore the dependence of our estimates on the GPDs we compare in Fig. IT2| the ratio R 
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Figure 12: The ratio R defined in ( ^31) for ^/s = 5 GeV and Q'^ = 5 GeV^. The curves 
correspond to the three models if^ = Hf)^ (sohd), H'^ = Hj^j^ + Hjj (dotted), and 
ff" = Hf)^ - Hi (dash-dotted). 



for the cases where H'^ is taken from the double distribution ansatz (|3^) , (pq ) alone, or 
as the sum of this and the D-term in (^). Due to a numerical accident the contributions 
from -f/^DD, and Hdd in ( p3D nearly cancel each other and produce a quite small inter- 
ference term. This result should be interpreted with care since, as we discussed, is 
obtained by extrapolating information from the usual parton distributions into the ERBL 
region, and our H'j-, is the result of a particular dynamical model. With a generic D-term 
of the same size one could also obtain a rather sizeable interference signal, as we see in 
Fig. |l^ when combining H'j^jj and H^) with the "wrong" sign for In accordance with 
our discussion at the end of Sect. ^Tl| we conclude from this exercise that the unpolarized 
interference term is highly sensitive to the behavior of the GPDs in the ERBL region, 
where our modeling is least reliable. 

Fig. |T3| shows R for the same three models of H'^, now as a function of r at fixed Q'^ 
and t, and thus for varying collision energy a/s. Notice that at r = 0.36 the minimum 
value of \t\ is equal to 0.2 GeV^ so that one is in coUinear kinematics, where the angle 
is undefined. As the numerator of R projects out the coefficient of a cos dependent 
term in the cross section, it must strictly vanish at that point. We remark that r = 0.36 
is still far from its maximum value Tmax = (1 + 2M/(5')~^ = 0.54, where the production 
threshold ^/s = Q' -\- M is reached for Q'^ = 5 GeV^. It is interesting to note that in 
TCS the total collision energy at threshold is large, whereas in DVCS or in inclusive 
deep inelastic scattering one scans the resonance mass region down to the proton mass as 
xb approaches its upper limit 1. While the straightforward application of leading-twist 
dominance seems dangerous in TCS at r close to Tmax, this might be an interesting regime 
to study parton-hadron duality. 

We finally wish to remark on TCS with a neutron target. In that case the BH process 
is suppressed due to the zero charge of the neutron. We can explicitly see this in the 
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Figure 13: The ratio R for the same models as in Fig. |12|, but as a function of r at 
Q'2 = 5 GeV^ and \t\ = 0.2 GeV^ 



approximation (P^D, where the term in brackets involving the (typically large) factor l/r^ 
goes with a combination of form factors that vanishes for t —>■ 0. We find that the TCS 
contribution to the cross section is indeed more important than for a proton target. With 
the kinematics in Fig. |1^ and the double distribution ansatz (|3^), (|36D it does however 
not amount to much more than 10% of the BH contribution. The unpolarized interference 
term on the other hand generates a tiny ratio R of barely 1%. This can be understood 
from (^), where the potentially large contribution from Tii is suppressed by the Dirac 
form factor. Tii, whose prefactor survives in the t —>■ limit, is penalized with a small 
factor rj and further suppressed by a near cancellation of the charge-weighted polarized u 
and d quark densities in the neutron. 

6 Summary and discussion 

Next to DVCS, timelike Compton scattering may be the theoretically cleanest process 
where generalized parton distributions can be accessed. To leading twist and at Born 
level, both processes involve in fact the very same integrals over combinations of GPDs. 
At the level of corrections and the departure from the large Q'^ limit, they will be 
different. A simultaneous description of both reactions may thus be a benchmark test of 
our understanding of the dynamics, both of the approximations employed in describing 
the parton level process and of the nonperturbative input. 

TCS can be measured in exclusive lepton pair production, either with quasi-real 
bremsstrahlung photons from incident leptons or with a dedicated real-photon beam. 
Unlike DVCS, timelike Compton scattering is always accompanied by a Bethe-Heitler 
contribution much bigger than itself. It offers however relatively simple access to the 
real part of the Compton amplitude via the angular distribution of the produced lepton 
pair. Appropriate angular observables allow a rather clean investigation of the detailed 
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structure of the Compton process. With circularly polarized incident photons one has 
access to the imaginary parts of the Compton amplitudes and thus to the timelike analog 
of what has been observed in the DVCS channel using the lepton spin asymmetry [Q. 

Using the quark handbag diagrams of Fig. § and simple models of the relevant GPDs 
we have estimated the cross section and angular asymmetries for lepton pair production 
in a kinematical setting typical of the HERMES regime. We find that the angular asym- 
metry carrying information on the Compton process ranges from about 5% to 15% within 
the variations of the GPD models we have explored. This rather wide range of predictions 
is generated by varying the GPDs in the ERBL region. It illustrates that the real part 
of the Compton amplitude is highly sensitive to the form of these distributions in the 
region where their physics is least known and most different from that of ordinary parton 
densities. Similar results have been obtained in recent studies of the lepton charge asym- 



metry in DVCS 0. Given the uncertainties in modeling detailed features of GPDs, 
the numbers we estimate here should hence be taken with due care. We also remark that 
substantial as corrections to the DVCS amplitude have been reported for our kinematics 
p5| , p6[ . In any case, whether the real part of TCS is observably large or not will already 
provide important information about the dynamics of the Compton process. 

A look at Eq. ( |55D and the plots in Figs. |^ and ^ reveals that in the region of r we 
considered, the imaginary part of the Compton amplitude is significantly larger than its 
real part. From Eqs. (^) and it then follows that in this kinematics the photon 
helicity asymmetry will be larger with our model GPDs than the unpolarized angular 
asymmetry we have investigated here. 

We have not attempted to give estimates for the regime of very small r, where DVCS 
has been observed at the HERA collider |^. In that case, the contribution from gluon 
GPDs at order as is expected to be too important to be neglected. As to the Born level 
quark contribution below r = 0.1, we find in our model that as r decreases both ReTii 
and — ImTii rise, as well as their ratio — ReTii/ImTii. 

Contrary to what one might expect, our estimates for a neutron target do not give a 
much enhanced TCS signal in HERMES kinematics, neither in the cross section nor in 
the angular distribution. This is due to an unfortunate combination of mostly kinematic 
prefactors in the formulae for an unpolarized target. On the other hand, we do not expect 
such a suppression for coherent scattering on a deuteron target, whose GPDs have recently 
been discussed [HH . 
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Appendix 



In this appendix we discuss to which extent the evaluation of generalized parton distribu- 
tions in the double distribution model ( |37D to ( p9D requires knowledge of the usual parton 
distributions down to x = 0. This is of practical importance since parton distributions 
are of course only constrained by data down to some finite value of x, below which one 
must rely on extrapolations. 

As is well-known, the terms in ( P7D going with 6{x — x' — rjy') and 6{x + x' — rjy') 
individually have non-integrable singularities at x' = 0, but their sum is finite. To make 
this explicit we consider the charge conjugation even combination h+{x,r]) = h{x,r]) — 
h{—x,T]), which can be written as 
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(45) 



with q-\-{x) = q{x)+q{x). Note that the corresponding singularities at x' = are integrable 
for the quark valence combination g(x') — q{x'), which is not needed in the evaluation of 
the Compton amplitude. The singularities for the polarized quark densities are integrable 
as well, and we can restrict our discussion to the most problematic case (^). Writing 
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(46) 



we see that the integrand in the first line of (^Sj) only involves x'g+(x'), whose singularity 
at x' = is integrable. 

The integrals in the second and third lines of (^) involve q+{x') down to values of 
order x' ^ rj — x and hence are potentially problematic for x ^ rj. To investigate them 
more closely, we decompose 
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(47) 



The first of the two terms in brackets leads again to the combination x'g+(x') and causes 
no trouble when the lower integration limit goes to zero as x rj. The second term 
does not provide a factor x' but a factor [r] — x) instead. Let us for the sake of argument 
assume that for small x the quark density behaves like 



with some A < 1. For x 



xq+[x) ~ X 
f] the integral involving the second term in 

— x\ f dx' qj^{x') ^rZ"^ —\ri — x\ 
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(48) 

then goes like 
(49) 
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and hence vanishes in the hmit where its evaluation requires knowledge of q+{x') down to 
x' = 0. 

To get a rough feeling for the integral in the first line of (^) and for the ones involving 
the first term in (|47D let us consider / dx' x'q+{x') with lower limit and upper limit of 
order rj. If we assume the power behavior (^HD, then the contribution from the interval 
x' e [0, e] to the total integral is of order {e/riY~^. For typical values of A this is about 
10% if e is one to two orders of magnitude smaller than t]. Of course q+{x) is unknown 
below some value of x, but unless its small-x behavior is much steeper than (^81) , the 
above estimate should not be altered significantly. 

Our discussion can be adapted to profile functions n{x', y') other than the one in (pQl). 
Provided that TT{x',y') is different iable in y' one can replace (^61) with 

7i(x', - n(x', = d27^(x', -) + 0(x'2). (50) 

where 821^ = dn {x' , y') / dy' . Decomposing 

x'jl -ri) r]~x 

= l-x 51 

77 V V 

we further see that (|47|) can be replaced with 
X — x'\ x'{l — rj) 
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d2TT{x', 1 - X') + 0(X'2) 



d27r(x',l-x'- '''^^ ^) )+0fa_^)2)^ (52) 



provided again that 7r(a;', y') is different iable in y' and that in addition 7r(x', 1 — x') = 0. 

Notice that while the GPD so obtained is finite at x = 77, its first derivative in x is 
in general not. For profile functions 7T{x',y') vanishing at x' + y' = 1 we readily obtain 
a representation of dh^/dx analogous to (^3|) with vr replaced by ?7^^927r. Provided that 
821^ is differentiable in y' our previous line of arguments goes through, except for the 
equivalents of (|47|) and (p2|) . There one will have an extra term T]~^d2'n'{x' ,1 — x') on 
the right-hand side, which for our profile function (|39D is nonzero. This means that the 
integrals corresponding to the second and third lines of (|45D will have an integrand going 
like q+{x') at x' — > 0, without a factor [rj — x) in front, dh^jdx thus behaves for x ^ rj like 
/ dx' q^{x') with lower limit of order \ri — x\, and with the small-x behavior (|48|) diverges 



hke 



-^h+{x,r]) 



''hr]-x\-\ (53) 
A 



Now, the principal value integral over x which gives ReTii according to (|^) effectively 
involves dh+/dx at x around rj, but the singularity ( ^3|) is integrable and gives a finite 
result for the amplitude. One can also insert the expression (|37|) into (§) and explicitly 
carry out the integrals over x and y' for the profile function (|39|) . The result for ReTii 
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has the form Jq dx' x'q^{x') p{x',ri), where p{x',ri) has a log(x') singularity at x' = 0. Up 
to this logarithm the small-x' behavior of the quark density thus enters Re Tii in the same 
way as according to our above discussion it enters h+{ri,ri) and hence Im?-^i. 

Let us explore how our arguments work at the quantitative level, restricting ourselves 
to u quarks for simplicity. In Fig. |n| we plot h^{x,ri) at rj = 0.11 as it is obtained from 
(|45|) when setting u^{x') = for x' below some cutoff e. We find good convergence as 
e approaches zero, the curves for e = 10~^ and 10~^ being hardly distinguishable. All 
predictions in this paper have been obtained with e = 10~^. Fig. |l^ shows Re?i" and 



ImTi]* as functions of r, calculated with the same cutoffs as in Fig. ^ when constructing 
the GPDs. Clearly the convergence is much slower for the real part, which can be traced 
back to large cancellations in the relevant integrals.^ This illustrates again the sensitivity 



of Re Til to small changes of the GPDs we have already encountered in Sect. 5.1 



In Fig. |T6| we show h^{x,r]) obtained from different parameterizations of the usual 
quark densities at jj? = 5GeV^, calculated with a lower cutoff e = 10~^ on x'. We 
compare the GRV 94 LO parameterization used in our predictions with three NLO dis- 
tributions in the MS scheme: GRV 94 NLO P7|, GRV 98 NLO [09|1, and MRSA' 



All input densities are clearly distinct for x below 0.01, but above x ~ 0.05 the three 
NLO parameterizations hardly differ among themselves. The corresponding curves for 
h'^{x,ri) at 7] = 0.11 are almost identical for the NLO parameterizations, in accordance 
with our arguments about the relevance of small x in the input densities. In Fig. O we 
make the same comparison for KeTi.^ and Im?i". In line with our previous findings, dif- 
ferences between the parameterizations are more prominent for ReTi", but they remain 
quite small. We thus do not confirm the results of pB, 1411], where substantially different 



GPDs and Compton amplitudes have been obtained from the GRV 98 NLO and MRSA' 
distributions with the same ansatz (|37|), ( p9| ) we have used here.|^ 

In summary, we have discussed in detail the calculation of GPDs with the double 
distribution ansatz (^7\). In the DGLAP region the calculation of h+{x,7]) involves the 
usual quark densities down to x' = {x — rj)/[l — ri), whereas in the ERBL region the ansatz 
literally does require knowledge of q+ix') down to x' = 0. The corresponding integrals 
involve however only the combination x'q+{x'), or they vanish like p9|). The contribution 
from at values of x' several orders of magnitude below rj should thus be moderate, 

unless one has an extremely steep rise of the quark density at small x'. We confirm this 



^°For the numerical evaluation of (||) we add and subtract H{x,r],t) at the points |a;| = r). To avoid 
cancellations as much as possible we only do this for |a;| < 2r] when rj < 1/2, writing 



PV / dx 
Jq 



X ~ 1] X + r] 



2r, 

= I dx 
'o 

■ 1 

dx 

'2i) 



1 



1 



X — rj x + rj 
1 1 



a; — 77 X + 7] 



h+ix,v)-h+ (77, 77)J + /i+ (?7, '/) In - 
h+{x,T]), 



where PV denotes the principal value prescription. 

^^Our results do not change qualitatively if we implement the double distribution ansatz 



H = 5 GeV but at ^ — 4 GeV as was done in 
under investigation with the authors. 



111. The origin of our discrepancies with 



not at 
56l |l| is 
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Figure 14: h'^{x,r]) for rj = 0.11, corresponding to r = 0.2, evaluated with a lower cutoff 
e = 10~" on x' in the integrals 
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Figure 15: ReTi" (left) and ImTi" (right) divided by ^F^{t), calculated from the GPDs 
in Fig. [14. 
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Figure 16: Left: Different parameterizations of u{x) + u{x) for fj? = h GeV^. Right: The 
corresponding generahzed distributions h^{x,ri) at r] = 0.11, obtained with our double 
distribution ansatz. 
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Figure 17: ReTi^ (left) and ImTi" (right) divided by ^F^{t), calculated from the same 
parameterizations of u{x) + u{x) as in Fig. ITB|. 



28 



numerically by evaluating the integrals (^) with different lower cutoffs on x' . For 77 of 
order 0.1 we find rather similar GPDs when implementing the ansatz ( P7D with different 
parameterizations of the quark densities, given that they show only mild discrepancies 
down to x' of order 10~^. 

The derivative of h+{x , rj) at x —>■ 1] becomes infinite with the profile function we used 
in our double distribution ansatz, but the corresponding singularity in x is integrable. 
The small-x behavior of q+{x) is found to affect the real and imaginary parts of the Born 
level Compton amplitude in a similar way, with the real part showing somewhat higher 
sensitivity. Again we confirmed this in our numerical study. 

We finally emphasize that our discussion of how relevant the usual quark densities at 
very small x are in the construction of GPDs refers to a particular model prescription. 
It is a different question to what extent one physically expects the behavior of q+{x) at 
X — *^ to be reflected in GPDs at finite rj. In the representation of GPDs as the overlap 



of wave functions for the incoming and outgoing hadron ||4^ one finds indeed that for 
X ^ 1] a. quark momentum in one of the two wave functions goes to zero. This is similar 
to, but not the same as the situation for a usual quark density at a; — 0, where both wave 
functions involve a quark with zero momentum. 
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